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Abstract 

Using a parabolic equation, we consider ray propagation in a waveguide with the sound speed 
profile that corresponds to the dynamics of a nonlinear oscillator. An analytical consideration of 
the dependence of the travel time on the initial conditions is presented. Using an exactly solvable 
model and the path integral representation of the travel time, we explain the step-like behavior of 
the travel time T as a function of the starting momentum pQ (related to the starting ray grazing 
angle xo by pq = tanxo)- A periodic perturbation of the waveguide along the range leads to wave 
and ray chaos. We explain an inhomogeneity of distribution of the chaotic ray travel times, which 
has obvious maxima. These maxima lead to the clustering of rays and each maximum relates to a 
ray identifier, i.e. to the number of ray semi-cycles along the ray path. 
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INTRODUCTION 



The ray tracing has demonstrated that early parts of the time-fronts, i.e. ray arrivals in 
time-depth plane, exhibit surprisingly regular structures simultaneously with sensitivity to 
the media inhomogeneities that give rise to ray chaos [1-5]. Even at very long ranges, up to a 
few thousand km, steep chaotic rays can form remarkably stable segments of the time-front, 
which are very close to their prototypes in the unperturbed waveguide. It turns out that 
each stable segment in the perturbed waveguide, like its counterpart in the unperturbed 
one, is formed by rays with equal identifiers, i.e. the number of ray semi-cycles. 

This phenomenon can be observed in the travel times of the so-called eigenrays, i.e. rays 
passing through a fixed observation point. In Ref. [6] (see also Refs. [7, 8]) it has been 
demonstrated that the travel times of chaotic eigenrays usually come in clusters with small 
time spreads centered at arrivals of unperturbed eigenrays. Although the rays that form 
a cluster have the same identifier, i.e. the same topology, this does not imply that their 
trajectories follow close paths. On the contrary, chaotic eigenrays contributing to the given 
cluster may significantly deviate from each other and from the unperturbed eigenray with 
the same identifier [3]. These results suggest that while the travel time of chaotic rays is 
a random function of starting parameters, it is much more predictable as a function of its 
identifier and the trajectory cndpoints, and it also relates to the dependence of the travel 
time T on the starting momentum po [9]. Interesting, even puzzling features, when observed 
numerically, reveal a step-like behavior of T" as a function of the initial momentum that 
forms the so-called "shelves" for ray propagation in a range-independent waveguide. In 
the case when rays propagate in a range-dependent waveguide so-called, ray chaos takes 
place, and an important characteristic is a ray travel time distribution, which exhibits fairly 
inhomogeneous features with obvious maxima [9]. 

In this paper we study analytically the dependence of the travel time on the initial 
conditions in the framework of an exact solvable model. We consider ray propagation in 
a waveguide with the sound speed profile which corresponds to the dynamics of a quartic 
oscillator. Therefore, studying this model, we explain the step-like behavior of the travel 
time T as a function of the starting momentum po (related to the starting ray grazing angle 
Xo by Pq = tanxo)- For the case when ray chaos takes place due to a range dependent 
perturbation, we also explain the inhomogeneity of the ray travel time distribution which 
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has obvious maxima. These maxima lead to the clustering of rays, and each maximum can 
be related to the corresponding identifier of the cluster of rays. 

The paper is organized as follows. In Sec. 2 we give a brief description of ray theory 
in the small-angle approximation. Explicit expressions for the Hamiltonian determining 
the ray motion and for the ray travel time are presented. An overview of the numerical 
results on the ray travel time, obtained in [9] , will be presented and explained in the present 
paper. An exact solution for a simplified speed profile corresponding to the quartic oscillator 
is considered in Sec. 3. An exact expression for the corresponding classical action as a 
function of the initial conditions is presented. A detailed analytical analysis of the step- 
like functional dependence of T on Pq is performed as well. In Sec. 4 the maxima of the 
distribution function for the ray travel time are found for the integrable quartic oscillator in 
the presence of a perturbation. This analysis corresponds to the semiclassical consideration 
of an amplitude of the wave function in the momentum representation. The conclusion is 
presented in Sec. 5. Some details of calculations related to the step-function are presented 
in Appendices A-C. 



RAY TRAVEL TIMES 



Parabolic equation approximation 

Consider a two-dimensional underwater acoustic waveguide with the sound speed c being 
a function of depth, and range, r. The sound wave field as a function of r, z, and time, 
may be represented as 

where the Fourier components u are governed by the Helmholtz equation (see for example [10, 

_ + _ + ,V„ = 0. (2) 

Here k = uj/cq is a wave number, while n = CQ/c{r,z), is the refractive index and Cq is a 
reference sound speed. For the 2D picture, grazing angles are defined as the ratio between 
the wave numbers kz and kr'. tanx = kz/K, where k — ■\/k1-\- k^. In the small-angle 
approximation, when sound waves propagate at small grazing angles with respect to the 
horizontal, i.e. kr ~ k, the Helmholtz equation may be approximated by the standard 
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parabolic equation [1, 10, 11]. Present u as 

u{r, z, uj) — v{r, z, u) e' 



(3) 



and substitute this expression into Eq. (2). Taking into account that u is a slowly-varying 
function of r and neglecting the second derivative of v with respect to r, we derive the 
parabolic equation 

This equation coincides formally with the time-dependent Schrodinger equation. In this 
case the partial derivative with respect to z is an analog of the momentum operator, i.e. 
p — —ik~^d/dz, while r plays the role of time and associates with the Planck constant. 
In underwater acoustics it is always possible to choose the reference sound speed Cq, such 
that |n — 1| <C 1, and replace 1 — by 2(1 — n) = 2(c(r, z) — cq)/cq. 

Since r is a time-like variable, the Hamiltonian system formally coincides with that de- 
scribing a mechanical particle oscillating in a time-dependent potential well U with the 
Hamiltonian 

i/ = ^ + C/(^), (5) 



where 



U{r,z)^ ■ (6) 

Co 



The dimension variable p is an analog to the mechanical momentum. It relates to the ray 
grazing angle X by p = tan^. The "potential" U in Eq. (6) represents a potential well 
whose parameters may vary with the range r. 

For the point source located at r = and z — zq we have 

where the sum goes over contributions from all rays connecting the source and the observa- 
tion point (r, z). Such rays are called the eigenrays. Here S{z, Zq, r) is the eikonal analog to 
classical action or the Hamilton principal function in mechanics - of the i^-th eigenray. This 
function is defined by the integral [12] 

J (pdz-Hdr) (8) 



over the ray trajectory from (0, zo) to (r, z). 
The amphtude A{z, zq, r) is given by [13] 



dzdzn 



\dz/dpoy 



(9) 



where C{uj) is a function determined by the time-dependence of the radiated signal, and the 
derivative dz/dpo is taken at the range r. 

Substitution of Eqs. (3) and (7) into Eq. (1) yields 



u{r, z, t) = / duj Ai,{z, Zo, r, a;) exp (iuj ( 1 Si,{z, zq, r) — t 

~f J V Vco Co 



(10) 



Each term in this sum represents a sound pulse coming to the observation point through a 
particular eigenray The quantity 



T ^- + -S{z,zo,r) 

Co Co 



(11) 



determines a delay of the pulse relative to an initially radiated signal and it is called the ray 
travel time. 



Numerical results: an overview of [9] 

Studying the general properties of ray travel times in acoustic waveguides is equivalent 
to studying the properties of the principal Hamiltonian function 5" of a mechanical particle 
oscillating in a potential well. Recently the properties of S have been numerically studied 
in [9]. Hereafter, we refer to this result as SVZ. The main numerical results important for 
the present analysis are shown in Figs. 1 and 2, which were taken from SVZ. Both figures 
present the travel time dependence on the starting momentum pq. Figure 1 demonstrates 
dependencies of the ray travel time T on the starting momentum po for two waveguides with 
different sound speed profiles, i.e. for two "potentials" U {z) : 

Ui{z) = az'^ + bz^, and U2{z) = --fz. (12) 

All six curves shown in Fig. 1 present the travel times at a range of 150 km, and each curve 
corresponds to a particular source depth. Even though the "potentials" Ui and U2 are quite 
different, both dependencies T(pq), shown in Figs. 1 and 2, have an important common 
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feature: each curve has "shelves" where its inchnations with respect to the horizontal are 
small. At intervals of starting momenta corresponding to the "shelves" , the ray travel time 
T is most stable (least sensitive) with respect to small variations in po- The same features 
occur for the so-called canonical sound speed profile or the Munk profile c{z) = cm{z), 
widely used in underwater acoustics to model wave transmission through a deep ocean 
[10, 14]. The dependencies T{pq) are shown in Fig. 2 which presents the ray travel times, T, 
as a function of the starting momentum, po, at range of 2500 km for a point source located 
at 1 km depth. A thick solid line graphs the T{po) for the regular ray in the aforementioned 
range-independent waveguide. Randomly scattered points correspond to chaotic rays in the 
presence of perturbation: c{z) = Cm{z) + 5c{z,r). The density of these points has maxima 
forming dark horizontal stripes, as is shown in Fig. 2. It should be pointed out that, while 
the background profile c{z) is realistic, the perturbation 5c{z, r) has been chosen in SVZ 
to present only a strongly idealized model of internal-wave-induced sound speed variations 
[15-17]. Nevertheless, this perturbation causes a chaotic ray behavior whose properties 
closely resemble that observed in more realistic numerical models [1, 2]. Ray travel times 
in the unperturbed waveguide {5c = 0) presented in Fig. 2 (thick solid line) have the same 
properties typical of range-independent waveguides, namely "shelves" , which are similar to 
those presented in Fig. 1. Moreover, for the chaotic rays, the stripes of the scattered points 
are located at travel times close to that of the "shelves" on the unperturbed T{po) curves. 
Note that the unperturbed "shelves" may be interpreted as parts of the T{po) curve with 
the highest density of points depicting unperturbed arrivals. It has been stated in SVZ that, 
under conditions of ray chaos, the positions of maxima of the density of ray travel times 
remain relatively stable. So, these figures express the main puzzling results of the generic 
features of "shelves" , and our main task in the present work is to explain them. 

A QUARTIC OSCILLATOR 

As was mentioned above, the general properties of the ray travel time T can be described 
by the action 5" of a mechanical particle in a potential well. Therefore, the generic features 
of "shelves" for an unperturbed ray can be explained in the framework of an analytical de- 
scription of the Hamiltonian principal function or the classical action S (8) for an integrable 
system with some potential U (6). Below we consider the oscillating dynamics of a particle 
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in the potential Ui of (12). 



Action 

As the momentum p = tanx is a dimensionless variable, it is convenient to consider 

the potential in the dimensionless variables as well. Namely we consider \plaz — > z and 
-\/2ar cijr, while hJ^ jc? = X. Therefore, the dynamical equation for a particle in the 
potential Ui (also called a quartic oscillator) in the new notation is 

z + u;'^z+ Xz^ = 0. (13) 

Following to SVZ we take a mass m — 1. We also use here, formally, the notation z = 
(Pz/dr^, i.e. the range r plays the same role as a formal time in Hamiltonian dynamics. 
This equation can be solved exactly. The solution is chosen in the form of the Jacobian 
elliptic function [18, 19] 

z{r) = Zcn{nr + (I)o,k), (14) 

where Z and 0o are an amplitude and an initial phase of oscillations correspondingly. The 
frequency of the nonlinear oscillations is 

n^^uj^ + xz^ (15) 

and the modulus of the elliptic functions is 

= X{Z/nf. (16) 

These values are obtained by the straightforward substitution of the solution (14) into (13). 
Following [19] we take the modulus k and the initial phase 4>o to be constants of integration. 
In this case, the solution (14) is 



z{r) = Zen {4>, k) 



2, ,2 



A(1-2k2) 



1/2 



cn 



+ 00, 1^ 



VI - 



(17) 



where k and 0o are associated with the initial coordinate Zq and momentum pq as 

Zq — z{r = 0) = Zca. (0o, >^)-i Po — z{r = 0) = —ZD,sn {(pg, K)dn (0o, (18) 

with sn and dn are also Jacobian elliptic functions. It is simple to see from (18) that k is 
the integral of motion related to the Hamiltonian 

K = ^/XH/ou^ (19) 



while the initial phase is 

(t)o^cn-\ujz/V2H]. (20) 

It also follows from (18) that for po > 0, the initial phase changes in the range 3A'(/t) < 
(f)Q < 4K{k) (or —K{k) < 4>q < 0), where K{k.) is the elliptic integral of the first kind. The 
modulus is restricted by < < 0.5, and the relations between the constants of integration 
and the initial conditions are expressed by the single-valued functions. 

Inserting (17) in (8), and using the integrals (312.02), (312.04), (361.02) of the Ref. [18] 
and the formulas for the elliptic integral of the second kind [19] 

E{(l)) - E{(f)') = E{(l) - (f)') - Khn (0)sn ((/)')sn (0 - 0'), 

we obtain the following expression for the action S 

S = Z^£;(Qr) + QV(1 - «;2)(2 - 3^)/3X + ^^^^^ x 

< sn (0o)sn (0)sn {fir) + ^[sn (0o)cn {(pojdn (0o) - sn (0)cn (0)dn (0)] ^ (21) 

where 

UJ 



fir + ^ 



VI- 2^2 ■ 

The following notations E{x) = E{x,k) and sn (x) = sn (x, k) (the same for cn,dn) are 
used. 

"Shelves" in the small n approximation 

The small k, approximation. The expression for the action S can be simplified. Since < 
0.5, one can use the small-AC-approximation for the elliptic integrals. Using the following 
definition for the elliptic integral [21] 

E[x, k) = E[x) = X — sn (x')dx' 

Jo 

and the approximation sn (x) ~ sin(x), we obtain approximately that 

E{x) fti X - xK^/2 - sm{2x) /4. (22) 

Inserting(22) in (21), and then combining the first two terms, we obtain, after doing small 
algebra, the following expression for the action 

^^^^ ^ " ^ ['^"^^'^^ " '^^^"^'^^ ' ^^^^ 



where the nonhnear frequency is now ft ^ 00(1 + k^). It also follows in this approximation 
that the relation (18) between the initial momentum po and the modulus k is simplified 

Po ~ gi^, (24) 

where g — —uj^/2/ Xsin 0o, and — 7r/2 < 0o < 0. The dependence of ray travel times 
on the initial momentum T{po) in SVZ coincides up to some constant multiplier with the 
dependence of the action on the modulus S{k,) in (23). 

"Shelves. " It follows that the action in the form (23) consists of two terms. The first one 
is the dominant (proportional to cur ^ 1) monotonic growing function in k. The second one 
is the small but fast oscillating term with a large frequency (proportional to a;r ^ 1). Such 
a combination of these two terms ensures the monotonic growth of the function in general, 
but at the same time the extrema equation dS/dn — has solutions. These solutions can 
be simply obtained, e.g. for (f)() = 0. The extremum points condition gives, in the same limit 
Lur ^ 1, the following solutions for k 

20 = 2Q{K,)r + 200 = ± arccos(2/3) + 27rm + 0{l/u;r) = 0^, (25) 

where m > ot/tt are integer numbers and 0(1 /uur) means neglected terms of the order of 
1/ujr < 1. 

Therefore, there are regions between extrema points (0~ , 0+ ) with the same number m 
where the derivatives are negative, dS/dn < 0. It follows that, in a range of An — A_ 
n/SLorn, the action S decreases by AS- (see Appendix A). These regions alternate with 
regions of growth, where dS/dn > 0. Between extremum points (0^, 0^+i) on the range of 
Ak, = A_)_ = 3A_ the action changes as 

AS+ = 9AS- (26) 

Therefore the growth of the action is stronger (by 9 times) than the decrease that leads to 
the step-like behavior of the action as a function of k. This step-like function (see Figs. 1 
and 3) has horizontal parts called "shelves" in [9]. 

An important feature of "shelves" is a large number of Fourier components in the Fourier 
transformation of the oscillating term in (23) (see Appendix B). It is shown in Appendix C 
that the average number of "harmonics" of the Fourier transformation is 

(Ds) «i u;r > 1 (27) 

One can see, in the insert of Fig. 3, a large number of the Fourier amplitudes. 
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TRAVEL TIME DISTRIBUTION FOR CHAOTIC RAYS 



In contrast to the regular dynamics, the arrival times of chaotic rays are not uniquely 
defined functions of the initial conditions, which is a simple result of the energy H (5) or the 
modulus K being no more the integrals of motion in the chaotic case. This means that many 
initial conditions can contribute to the same arrival time (as it is seen in Fig. 2). Wave 
dynamics leads to the wave superposition with different probabilities for different arrival 
times. Obvious maxima of the travel times distribution are seen in Fig. 2. To explain this 
phenomenon, we will use the analytical solution for the unperturbed ray dynamics, while 
chaotic dynamics is modeled by a randomization of the initial phase 0o or by a variety of 
different sources with random phases 0o uniformly distributed in the interval (— 7r/2, 0). 

An integrable case. The probability of finding a particle with the range r and the depth 
z is defined by a solution of the parabolic equation (4) with an amplitude (9). Therefore 
these amplitudes define the probability distribution for different S{k) with the same fixed r 
by (9) 

\A{r,z)\'cx\dz/dpo\-\ (28) 

Taking into account the solution z{r) in (17) and the relation Pq{k.) in (24), we obtain in 
the small-K-approximation 

dz/dpo = [dn/dpo) ■ {dz/dn) + (90o/9po) ■ {dz/d4>o) 

[AourK^ cos 00 sin — 2 cos(0 + 0o)]/ sin(20o) (29) 

In the limit ur ^ 1, the main contribution to the derivative (29) is due to the linear term 
(f) ~ cur. Therefore the evaluation of the probability for the asymptotically large times is 

1 



\A\ 



sm[(po) 



sin^ 



(30) 



The maxima of this probability correspond to zeroes of the denominator, which can be found 
from the following equation 

(p{K = Kn) = (po + ^{Kn)r = (l)o + u}r{l + kI) = Trn, n = 0,1,2,... (31) 

For the fixed cur and (po, the solutions of (31) /t = k„ determine the actions Sn — S{Kn) 
where the maxima of probability take place for the integrable case. 
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Ray chaos. For the chaotic case the energy H or the modulus k, are no longer the integrals 
of motion. Therefore the rays with different initial conditions k, 0o can contribute to the 
same arrival time S with different probabilities. In our phenomenological approach, it is 

convenient, as was mentioned above, to model the chaotic dynamics by a variety of initial 
conditions with random phases 0o uniformly distributed in the interval (— 7r/2, 0). Therefore 
the averaged probability is a superposition of probabilities with all initial phases. It reads 

(UP) = - / \A\^d(Po = (sm(Qr)/8 - (I/tt) cosfl^r) In [- tanfl^r)]) , (32) 

where the signs ± are due to the modulus function in (30) and (32), and (+) sign stands 
for — 7r/2 < (Qr, mod 27r) < 0, while (— ) sign is taken for 7r/2 < (Qr, mod 27r) < tt. The 
maxima of the mean probability are 

^{i^n)r — (^r{l + kI) — TTO, (33) 

that coincides with (31) for 0o = 0. It follows from (31) and (32) that rays with different k 
are clustered by the index n numbering the maxima. For all values of 0o one always finds a 
value of K which corresponds to the maxima conditions with the same n. It also follows that 
all other values of k which do not correspond to the maxima conditions "carry" the same 
index n if their action S is close to the maximum value 5'(«;„). This phenomenon of the ray 
clustering can be a possible explanation of the ID number for rays [9] . 

CONCLUSION. 

It should be admitted that in the framework of this simple analysis of the solution of the 
quartic oscillator, we are able to describe fairly well the step-like behaviour of the arrival 
paths as a function of the initial momentum. This step-like behaviour is known as "shelves" 
[9]. 

For the chaotic behaviour of rays, we constructed a phenomenological model and pre- 
sented only qualitative explanations of the nonuniform distribution of the arrival paths as 
a function of the initial momentum. The maxima of this distribution are explained in the 
framework of the integrable model. Such a kind of consideration corresponds to a so-called 
linear response approximation. 

This work was supported by the U.S. Navy Grant N00014-97- 1-0426. We thank A. Pitt 
for her help in preparing this paper. 
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Appendix A 

The extremum points condition gives in the same hmit cur ^ 1 the following solutions 
for K 

2(f) = 2Q{K)r + 200 = ± arccos(2/3) + 27rm = 0^, (A. 1) 

where m > ujr/'K are integers numbers. The phases 0+ stand for the minima of S with 

d''S{(t)^/dK^ = S';^^^ = sVbJ'nS-^ > 0, (A. 2) 

while 0^ define the maxima of the action, 

d''S{(l)-/dK' = S'l^_ = -8v^a;^KV'. (A. 3) 

It is simple to see that the regions on k between any adjoint extrema are very small. Indeed, 
the width of the region where the action decreases, A_ is determined from (25) 

0+ - 0- = 2cjr[{K + A_)2 - K^] 7r/2, 

where we took approximately that arccos(2/3) k, ^tt/A. From where we obtain that 

A_ 7r/8ujrK. (A. 4) 

Analogously, from 0~_,.i ~ 0m ~ 37r/2 we obtain that the width of the region where S 
increases is 

A+ = 3A_. (A. 5) 

Since A± ^ 1, we can define both a growth AS'+ and a decrease AS"- of the action in corre- 
sponding regions between adjoined extremal points in the linear approximation. Expanding 
the first derivative dS/dn near every extremal point, we obtain for /S.S^ 

/■A-/2 I-A-/2 

A5_ = / S':^^_xdx + / S':^^^{-x)dx. 
Jo Jo 

Inserting (A. 2) and (A. 4) in the integration, we obtain that 

AS- = -Ti^V5uj^Kyi6. (A. 6) 

Carrying out the same for AS+ we obtain 

A5+ = 9|A5_|. (A. 7) 
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Appendix B 



Let us rewrite the oscillating term in the form 

sin(20) — sin(2a;r + 20o) cos(a;r/«^) + cos(2a;r + 20o) sin(a;r/«^). (B. 1) 

For simplicity we consider k e [0, 1] by rescaling 2/t^ that docs not lead to any errors 

in the analysis. Since the region of definition of sin(20) is restricted by this segment, it 
is not difficult to show that the coefficients of the Fourier transformation /'^(s), /'^(s) are 
determined by the Fresnel integrals C{s),S{s) [20, 21]: 

/^(s) = / dKsm{2(f))cos{27rsK), /^(s) = [ o?«;sin(20) sin(27rsK). (B. 2) 

Jo Jo 

Carrying out the variable change x — ujtk and considering that or ^ 1 we take the upper 

limit to oo. Then we have for (B. 2) the following four integrals which determine the 

coefficients f'^{s),f^{s) 

pi 1 s 

/ dKsm(u!rK^)sm(sK) ^ —= / (ixsin(2;^/a;r) sin( — x) 
Jo V'^r Jo ujr 

and it gives 



TT I f S \ 

cos C I — ^= + sin S 



2ujr \ Aujr \2^/ujr J Auor \2^/LJr ^ 

Analogously we obtain for the rest of integrals 

dKsm(ijjrK^) cos( sk) ~ \ I < cos sin - — /- , 

Hcur { Acur Acur J 



lo 



A 



dKcos{urK ) cos(sk) ^ J - — < cos h sin 



V 8a;r [ Aujr Aujr 

I dKcos{ujrK^)sin[sK) 
Jo 



TT I f S \ 

sin - — C — — cos - — S 



2ujr \ Aur \2^/oJr J Aur \2^/uJr ^ 

Keeping AsAx > const, we obtain that there are of the order of or ^ 1 components with 
the amplitudes ~ 1/ ^/uJr contributed to the Fourier transformation. 
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Appendix C 

The oscillating part of the action S has a complete oscillation between points (0^, (l^m+i) 
that corresponds to the range on k or a quasi-period equaled to 

D« = A+ + A_ = 4A_. (C. 1) 

Hence, taking into account (A. 4), we obtain that the number of harmonics in the Fourier 
transformation is 

D, = 2n/D^ = iujrK. (C. 2) 

Since < k.^ < 0.5, the averaging of (C. 2) gives 

{D,} = ojr. (C. 3) 

It should be stressed that this estimate is approximate and gives only the order of Dg. 
The exact theorem on the uncertainty conditions (see for example [22]) ensures only that 
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Figure captions 

Fig. 1. The ray travel time as a function of starting momentum for two waveguides with 
the sound speed profiles Ci{z) — Cqi + az'^ + bz^ (curves a, b, and c), and C2{z) — C02 — "yz 
(curves d, e, and /). Parameters: Cqi = 1.49 km s~^, a — 1. km~^s~^, b — 1. km~^ s~^, 
C02 = 1.4465 km s~^, 7 = 0.0435 s~^. It has been assumed that the waveguide with ci{z) has 
no boundaries, while C2{z) has a reflecting surface at 2; = 0. The travel time at each curve 
is estimated from the arrival of the ray with po — 0. Different curves present rays escaping 
point sources located at depths: km (a), 0.5 km (6), 1 km (c), km (d), 1 km (e), and 2 
km (/). [from Ref. [9]] 

Fig. 2. The ray travel time versus starting momentum in the unperturbed (thick solid lines) 
and perturbed (points) waveguides at the range of 4500 km and for the point source set at 
a depth of 2.5 km. [from Ref. [9]] 

Fig. 3. The ray travel time (action S versus the modulus k. for Eq. (23), where 0o — 
— 7r/4, u; — 1, r — 355.4, A = 1.2. The insert is the amphtudes f{s) vs s of the discrete 
Fourier transformation for the oscillating part of the action S (B. 2). 
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